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We investigate the linear ac-conductance for tunneling through an arbitrary interacting quantum 
dot in the presence of a finite dc-bias. In analogy to the well-known Meir-Wingreen formula for the 
dc case, we are able to derive a general formula for the ac-conductance. It can be expressed entirely 
in terms of local correlations on the quantum dot, in the form of a Keldysh block diagram with four 
external legs. We illustrate the use of this formula as a starting point for diagrammatic calculations 
by considering the ac-conductance of the noninteracting resonant level model and deriving the result 
for the lowest order of electron-phonon coupling. We show how known results are recovered in the 
appropriate limits. 
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Introduction.- Recent experimental progress enables 
the study of the dynamics of electronic processes in 
mesoscopic structures. In a seminal experiment, Gabelli 
et al. [1^ measured the ac-conductance of a mesoscopic 
RC-circuit. Measuring the in and out of phase compo- 
nents of the ac-conductance at GHz frequency, they con- 
firmed theoretical predictions for the coherent transport 
of noninteracting electrons 13 Sj • General interest in 
time- or frequency-dependent properties of mesoscopic 
systems is driven by efforts to measure and manipulate 
quantum coherent systems as quickly as possible, for in- 
stance, for applications in quantum information trans- 
mission and processing. An important building block for 
such experiments may be provided by a coherent single- 
electron source recently realized Furthermore, time- 
or frequency-dependent measurements can yield addi- 
tional information about internal time and energy scales 
that are not accessible via the time-averaged dc-current. 

A cornerstone for many recent theoretical considera- 
tions of the time-independent dc-current was laid in the 
1992 paper of Meir and Wingreen [6 , who derived a 
Landauer-type formula, giving the current through a cen- 
tral interacting region in terms of local Keldysh Green's 
functions of the central region (and the tunnel-coupling 
and Fermi-functions of the leads). In this paper we ex- 
tend this approach to the finite frequency case and de- 
rive an expression for the current response to a small 
ac-voltage excitation, which is applied across the cen- 
tral interacting region in addition to a static dc-bias, 
i.e., the linear ac-conductance in a nonequilibrium sit- 
uation. As in the Meir-Wingreen approach we derive a 
formula for the current in terms of 'local' objects, but 
in our case these are no longer the usual Green's func- 
tions, i.e., expectation values of two-operator objects, but 
contain four electronic operators. The actual evaluation 
of these central objects may in principle be performed 
employing various techniques and approximations. The 
approach lends itself, however, to a combination with a 
diagrammatic calculation perturbative in the interaction 
strength , as we demonstrate in an examplary calcula- 



tion for electron-phonon coupling. 

Before turning to the general approach, we note some 
previous works that have considered special cases. For 
vanishing electron-phonon coupling our approach repro- 
duces calculations based on the scattering matrix method 
[21 in E]. Those calculations rely on a noninteracting 
single-particle picture, but can be amended to capture 
screening by incorporating the effects of induced charges 
on the internal potentials in a self-consistent manner. 
The current response to a fixed internal potential as cal- 
culated by the scattering method (or by our approach) 
is then used as an input for these calculations. Another 
route to finite frequency calculations starts from an ex- 
pression of the current in Green's function formulation 
for completely general time dependence of the parame- 
ters [9]. Progress can then again be made for the non- 
interacting case [9] , or employing special approximations 
to the interacting self-energy and Green's function |10j . 
Kubo relations naturally link the linear ac-conductance 
to current-current correlators. Indeed, as pointed out 
in [11^ its real part is related to the asymmetric part 
of the frequency dependent current-noise for a static dc- 
bias. Derivations of a general Meir-Wingreen type for- 
mula for the noise, however, have been restricted to the 
symmetrized noise and, furthermore, consider the nonin- 
teracting case |12j or focus on the zero frequency limit 
[13l[14]. 

System.- We consider transport through a central re- 
gion, for instance a quantum dot, between a left and 
right lead. While electrons in the leads can be consid- 
ered as noninteracting, within the central region there 
might be electron-electron or (as considered in the spe- 
cific calculations in this paper) electron-phonon interac- 
tion. The Hamiltonian of the system can then be writ- 
ten a,S H = i/ccntial + Y^a=L/RiHres,a + Hr.a) ■ Lead 

Hamiltonian -ffros,a = X)fe ^kacl^Cka and tunnel Hamilto- 
nian Hx.a = J2k n ^ka,ncl.^dn + c.c. are of the standard 

form and the central part, i?contrai = J2n^ndi_dn + i?int 
may include electron-electron or electron-phonon cou- 
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pling within Hint- We want to consider the nonequilib- 
rium situation, where on top of a finite dc-bias a small 
ac-excitation voltage of frequency Wac is applied to 
lead p. The resulting finite frequency current 1^ fiowing 
from lead a to the central region is given by the linear 
ac-conductance, laiwac) — ^Q/3(^^ac) ^^(wac) • Consider- 
ing linear response with respect to the ac-excitation, we 
can use Kubo formalism, but due to the finite dc-bias 
the expectation values in the Kubo-formalism have to be 
taken for the nonequilibrium finite-bias state. The stan- 
dard Kubo approach directly yields the linear conduc- 
tance Gai3 in terms of a (retarded) current-current corre- 
lator, fJc</3(Wac) = [Kai3{uJi,c) " Kap{0)] /(jt^ac) , whcrC 



Lit), ipiO) 



(1) 



The current operators are of similar form as the tunnel 
Hamiltonian, 

(— ej n n 

= y^XVka^nclgdn - Vn,kadn^ka) -(2) 



The correlator Kaf3{t,t') — K^pit — i',0) has the same 
structure as a retarded Green's function, but with a 
pair of creation and destruction operators for each of 
the current operators at the two time variables, see the 
schematic diagrammatic representation in Fig. [l] Each 
current operator appears as an external vertex with one 
lead-electron line (dotted line) and one dot-electron line 
(double line), one of them entering, one leaving the vertex 
[T5] . The two possibilities (for each vertex) carry signs 
according to Eq. In general, in between the four 

legs starting from the two vertices there will be a compli- 
cated diagram involving many interaction lines, as indi- 
cated by the shaded box in Fig. [T] The simplest type of 
diagram, however, is encountered in the non-interacting 
case. Then all operators are directly contracted, resulting 
in two Green's function running antiparallel (one starts at 
t going to i', the other goes from t' to t). Langreth's the- 
orem then yields so-called analytical continuation rules, 
giving the proper combination of nonequilibrium Green's 
function for the different Keldysh components of the 
product [7]. For instance, the object with two antipar- 
allel Green's functions C{t,t') = A{t,t')B{t' ,t) has the 
retarded component of C^{t,t') = A<{t,t')B^{t' ,t) + 
A^{t,t')B<{t,t') = liA^B"^ + A^B'^). These rules 
correspond to applying the standard rules for Keldysh 
indices at the outer vertices. (We use the following ma- 
trix representation of the nonequilibrium Green's func- 
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where the Keldysh 



Green's function component G^^ vanishes by construc- 
tion.) Turning to the interacting case, we use these same 
rules for the external vertices of the full correlator. We 
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Figure 1: Correlator needed for the ac conductance. One lead- 
electron line (dotted line) and one dot-electron line (double 
line) are connected to each outer vertex. The two contribu- 
tions in the current vertex, Eq. ([2|, corresponding to an in- 
or outgoing dot-electron line (and vice versa for the lead elec- 
tron) yield different signs for the various configurations. The 
four diagrams at the right show the possible combinations of 
Keldysh indices at the inner legs, with the lower two van- 
ishing. A single solid line stands for either dot-electron or 
lead-electron lines. 



find four possible combination of Keldysh indices on the 
inner legs of the correlator loop, see Fig. [T] Note that 
two of these combinations would include the vanishing 
Green's function component G^^ without further inter- 
action lines added to the simple loop diagram. In that 
case, the Langreth result is directly recovered. In the 
following we will argue, that the two latter index com- 
binations vanish also after including all interaction lines 
and we are left with only the upper two combinations in 
Fig.[l] 

To start with, we consider only one additional inter- 
action line. Then it is easy to see that this statement 
holds, simply by constructing all diagrams according to 
the indexing rules at the vertices. Vanishing diagrams 
are then found for two reasons. First, the application 
of the Keldysh indexing rules at the vertices may lead 
to the appearance of a line corresponding to the vanish- 
ing Green's function component G^^ . Secondly, diagrams 
may also vanish due to the fact that the appearance of 
an advanced or retarded Green's function between two 
vertices implies a fixed order of the time variables of the 
two interconnected vertices. Such a time order can lead 
to a logical contradiction, if, for instance, the Green's 
functions implied that ti <t2 < and ^3 <ti. 

To exploit this property systematically for all higher 
order diagrams, we introduce the notion of a time-arrow 
line, which appears in diagrams indicating a fixed time- 
ordering of the vertices along this line, see Fig. [2] Such 
a line runs parallel to G^'^ — G^ and antiparallel to 
G^^ — G^. Analysing the 12 different possible in- 
dex combinations at a vertex (we consider the electron- 
phonon case here, but the argument is general and relies 
on the Keldysh indexing rules at the vertices only) we 
find that if a time arrow-line runs into a vertex, there 
will always also be such a line leaving the vertex. In 
consequence a time-arrow line will never stop at an in- 
ternal vertex, but will run through a diagram until it 
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Figure 2: Time-arrow lines, running parallel to G^^ or an- 
tiparallel to G^^ , do not stop internally, but progress through 
any possible internal vertex. They may hit themselves, form- 
ing a closed loop, which indicates a logical contradiction in 
the time ordering, or end at the external vertices. Additional 
time-arrow lines (cf. the shaded box) may also exist and en- 
force vanishing of a diagram if they form a loop. 



hits an outer vertex (where it may end or go on) or un- 
til it hits itself and forms a closed loop. Such a closed 
loop will imply a logical contradiction within the time 
ordering as discussed above, and the corresponding dia- 
gram vanishes. Consider now a combination of Keldysh 
indices at the outer vertices as in the diagram shown in 
Fig. [2] Due to the vanishing of G^^ we can infer that two 
time-arrow lines start at the left external vertex, while at 
the right vertex we can have a time-arrow line starting 
or crossing through the vertex. As the time-arrow lines 
do not stop internally and can not stop at either of the 
two external vertices in this specific diagram, they have 
to form a closed loop. Analysing the external vertices in 
this manner, we can thus state that some index combina- 
tions will necessarily result in time-arrow loops and are, 
hence, forbidden. For the case of interest here, the two 
lower diagrams of Fig. [T] are of this type and vanish. 

Separating out the leads.- Each of the two external 
vertices is connected to a lead-electron line. Moving away 
from the vertices along these lines, we will encounter a 



tunneling event changing the lead-electron line to a dot- 
electron line [16]. Note, that there can be no interaction 
line in between the external vertex and the first tunnel- 
ing, as interaction is only present within the dot region. 
The Green's function corresponding to the considered 
line is now split into a product of a lead-lead Green's 
function, a tunneling matrix element and a dot-electron 
Green's function. The Keldysh indexing of those Green's 
functions follows the rules of matrix multiplication in 
Keldysh space, see Fig. [3] This split-off procedure has to 
be performed for the four possible orientations of dot and 
lead lines at the external vertices (taking account of the 
proper signs according to Eq. ([2|). We proceed by evalu- 
ating the sums over lead states contained in the outer cur- 
rent vertices and in the inner tunneling vertices marked 
by crosses in Fig.|3] For that, we change from time to fre- 
quency/energy domain. We introduce the frequency E, 
which encircles the loop, and the energy E' exchanged 
within the central object between upper and lower line 
(see Fig. |3|. The summations are then executed under 
the standard assumption of the wide-band limit to yield 

where ga^^^^ are the bare lead Green's functions, /" 
the Fermi-distribution function of lead a and is the 
couphng strength defined as F^ = 27r^j. jVap. For ease 
of notation we consider a single level in the central region 
here and in the following. The simple extension to the 
general case is discussed below. 

Separating out the leads in that manner we are left 
with a central object H with four legs of dot-electron 
lines, which includes all possible further interaction lines 
in its central block, see Fig. [3] Starting from the two 
remaining diagrams of Fig. [l] we can group together all 
contributions with the same indices on the central ob- 
ject n. We find that all but three contributions cancel, 
and the final result for the frequency-dependent retarded 
current-current correlator can be written as 
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Figure 3: Each correlator contains four combinations of in- 
/out-tunneling at the outer vertices. The dot-lead Green's 
functions are split into lead-lead and dot-dot Green's func- 
tions, so that the correlator can be expressed in terms of a 
central object 11 containing dot operators only, represented 
by gray boxes (with Keldysh indices at the legs indicated). 
For example, illi occurs in the first and last line, and these 
diagrams contribute to the term displayed in Eq. ([3|. 
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(3) 



X [?n? + (1 - 2/|+^,+^^j2n? - (1 - 2/|+^,)}n?] 

plus extra terms for the diagonal conductances, stem- 
ming from a direct contraction of lead operators of 
the outer vertices, indirect ^ _ e!r„ / |f [/|+^^^G^ -f 

Ie-u^G^] , so that = + Sc^pK^r' • 

This constitutes a general expression for the ac- 
conductance in terms of a central 'dot-operator only' ob- 
ject, which is the equivalent to the dot Green's functions 
in the Meir-Wingreen result. This central object H is 
a Keldysh-type diagram with four external legs. Two 
dot-electron lines enter the diagram and two leave it. 
We arrange the diagram so that the lines entering are 
placed at the upper left and lower right corner. Due 
to energy conservation the energies carried by the four 
external lines can be parametrized by three parameters 
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E,E', and Wac as indicated in Fig. [s] The four external 
legs carry Keldysh indices as indicated on 11. The central 
object n then contains all possible diagrams within these 
specifications, fully including tunneling and interaction. 

If we wish to consider several levels within the central 
dot-region, each dot-line carries extra level indices, in 
general difi'ering at start and end (due to tunneling to 
the leads and back into a different level). The couplings 
Ta become matrices in these dot indices and the trace is 
taken over the whole expression. For instance, the first 
term of Eq. ([3| would change into Ta^ijT/^^ki fi^ii with 
summation over all level indices i,j,k,l. In the following, 
we illustrate the application of our formula for two simple 
examples. 

Resonant level model- Without interaction, there are 
no phononic lines connecting the upper and lower line of 
n and energy exchange between these lines is not possi- 
ble. The superfluous E' integral is then cancelled by a 
(^-function guaranteeing energy conservation. The terms 
in Eq. ^ with and }nj vanish (since they would 
contain G^^ ~ 0), while the fHf contribution and the 
direct part immediately yield 
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The nonintcracting dot Green's functions including the 
tunneling self-energy are G^^^'^{E) = [E — e ± i{TL + 
Tfj)/2]^^ . We thus recover the result that can also be 
derived using methods specialized to the non-interacting 
case, viz. the fully time-dependent Green's function ap- 
proach of Ref . [9j or scattering matrix theory ^ . 

Second-order in e-ph coupling.- We now consider cou- 
pling to phononic modes of frequency lu,^, i.e. we take 
Hint = EK,[5K(aK + o.l)d^d-\- huj^ala^] . We find that 
in second order in the e-ph coupling g^, i.e., with one 
phonon line included in the diagram, there are three dif- 
ferent type of contributions. Firstly , we can replace one 
of the bare dot-electron Green's function by one dressed 
with a phonon line, see Fig. |4] For the evaluation of 
this type of diagrams, we accordingly replace bare by 
dressed Green's function in the zeroth-order result and 



expand to second order, e.g., G 



A,T 



G 



A, T^A, P^qA, T 
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where the self-energies YP^ are the lowest order Fock- and 
Hartree-terms, which depend on the bare phonon prop- 
agator Dk^^^^ . Secondly, there are contributions with 



one crossing phonon line, which yield 

o2 
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where the three terms in the square brackets correspond 
to the three different indexing variations of 11 in Eq. (13 
The contribution from the third type of diagram in Fig. 
vanishes in the dc-limit in linear response and will not be 
discussed further here. 






Figure 4: Three types of second-order diagrams. Diagrams 
with one Green's function dressed with a phonon line (left), 
diagrams with one phonon line crossing the loop (middle), 
and Hartree-type diagrams of two electronic loops coupled by 
a phononic line (right). 

The equivalence of the Meir-Wingreen result to the dc- 
limit of our result is most easily seen for the special case 
of linear response. Then the explicit difference of Fermi 
functions in the last two equations results in a deriva- 
tive, while all other Fermi and Green's functions can be 
taken at equilibrium and vanishing oj^c ■ Contributions 
from diagrams with one crossing phonon line can then be 
written 
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Expanding the Green's function in the Meir-Wingrcen 
result to second order in the electron-phonon coupling 
we finally recover the proper dc-limit. A slightly more 
involved calculation confirms the equivalence for finite 
dc-bias. 

Conclusions.- We derived a formula for the linear ac- 
conductance, Eq. ([s]), which is applicable for arbitrary 
interactions and for a finite applied dc-bias. The formula 
expresses the ac-conductance in terms of a central object 
n, which can be calculated using Keldysh diagrammatic 
methods. We illustrated the application of our formula 
by explicit calculations for the resonant level model and 
the case of electron-phonon coupling. 
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